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We report the current status of RBCK calculations on nucleon structure with both quenched and unquenched 
lattice QCD. The combination of domain wall fermions and DBW2 gauge action works well for isovector vector 
and axial charges, and moments of structure functions (x) q , di, and {l)sg- 
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Domain wall fermions (DWF) 1. have been 
successfully used by the RBC Collaboration for 
several applications In the baryon sector we 
had reported that DWF with Wilson gauge ac- 
tion in the quenched approximation can repro- 
duce the first negative-parity excited state of nu- 
cleon |3] which defied theoretical attempts for a 
long time. Here we report the current status of 
our calculations on the nucleon electroweak ma- 
trix elements, namely the isovector vector and ax- 
ial charges and some moments of the structure 
functions, in both quenched and full-QCD, using 
DWF and the DBW2 gauge action gj. 

Four form factors appear in neutron (3 decay: 
the vector and induced tensor form factors from 
the vector current, and the axial- vector and in- 
duced pseudo-scalar form factors from the axial- 
vector current. In the forward limit, the 4- 
momentum transfer should be small because the 
mass difference of the neutron and proton is only 
about 1.3 MeV. This makes the limit q 2 -> 0, 
where the vector and axial-vector form factors 
dominate, a good approximation. Their values in 
this limit are called the vector and axial charges 
of the nucleon: g v and g A . 

Lattice numerical calculation of these charges 
is relatively simple: prepare appropriate nucleon 
source and sink at rest, insert an appropriate cur- 
rent operator in between, and identify a plateau 
which gives a lattice estimate. There already 
had been several preceding works using Wilson 
or improved Wilson fermions p]. However since 
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Figure 1 . Quark mass dependence of axial to vec- 
tor charge ratio, g A /g v , quenched (circle) and full 
QCD (diamond). Full-QCD results, albeit with 
small statistics, are consistent with quenched re- 
sults: No mass dependence is seen. 



these fermions schemes do not respect the chi- 
ral symmetry, the lattice renormalizations of the 
relevant vector and axial currents are not the 
same. Worse, unwanted lattice artifact may re- 
sult in unphysical mixing of chirally distinct op- 
erators. The combination of DWF and DBW2 
preserves chiral symmetry well enough so that we 
are free of these difficulties: especially the rela- 
tion Z A — Z v is easily and well maintained, up 
to 0(a 2 ), or a few % at a « 0.15 fm. And in 
contrast to those preceding studies, we found a 
significant finite volume effect between the re- 
sults calculated on lattices with (1.2 fm) 3 and 
(2.4 fm) 3 volumes 0: about 20% increase was 
observed. On the large volume we find gA = 
1.212 ± 0.027(stat) ± 0.024(norm). The quoted 
systematic error is the dominant known one, cor- 
responding to current renormalization. This the- 
oretical first principles calculation, which does 
not yet include isospin breaking effects, yields a 
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Figure 2. The bare momentum fraction in 
quenched (circle) and full-QCD (diamond). 
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Figure 3. The bare (Vjsu-Sd matrix element in 
the quenched (circle) and full (diamond) QCD. 



value of gA only a little bit below the experimen- 
tal one, 1.2670 ± 0.0030 |7j. The dependence on 
the quark mass on the larger volume is very mild, 
while it is steep on the smaller. 

First few moments of nucleon structure func- 
tions are as easy as these charges, in that they 
do not require finite momentum states. The 
structure functions such as F\{x, Q 2 ), F2(x,Q 2 ), 
gi(x,Q 2 ), and g2(x,Q 2 ) are measured in lepton 
deep inelastic scattering, and h\(x,Q 2 ) in the 
RHIC Spin experiment. Their moments are fac- 
torized in the operator product expansion into 
Wilson coefficients and non-perturbative quan- 
tities such as (x n ) q , (x n )Aq, d„ or (i)sq- The 
latter are nucleon forward matrix elements of 
certain local operators, such as jD, 757Z?, and 
75CT-D, and hence are accessible by lattice QCD 
numerical methods similar to form factor cal- 
culations. There are calculations by QCDSF, 
QCDSF/UKQCD, and LHPC-SESAM collabo- 
rations with (improved) Wilson fermions [H]. 
In Lattice 2002 we presented some preliminary 
quenched results from the present DWF/DBW2 
calculations |Hj. The DWF/DBW2 combination 
indeed works better: good chiral behavior, ease of 
non-perturbative renormalization, and better ap- 
proach to the continuum were seen. We present 
higher statistics quenched results here, and some 
preliminary full-QCD ones. We still lack non- 



perturbative renormalization for some observ- 
ables. 

In the on-going full-QCD calculations, a 16 3 x 
32 x 12 lattice is used with the DBW2 gauge cou- 
pling set at (3=0.8 and the domain- wall height 
at M5=1.8. Dynamical DWF quarks with mass 
0.02, 0.03, and 0.04 (about strange mass) |5] in 
the lattice units are used. The nucleon structure 
calculations reported here are done only at the 
dynamical quark mass, and no partially quenched 
study has been made. About 900-200 conjugate- 
gradient iterations are necessary to achieve 10 -8 
accuracy for the Dirac equation solver. Hybrid 
Monte Carlo acceptance of about 70-80 % is ob- 
served. About 50 configurations are used at each 
dynamical quark mass value. The residual mass 
of m ros = 0.00137, p meson mass m p = 0.43(2) 
which translates into the lattice cutoff estimate 
of a -1 ~ 1.8 GeV, and the nucleon mass = 
0.56(4) or -1.01(7) GeV are obtained. Quality of 
plateaus in the three-point functions is similar to 
the quenched case at similar statistics. 

In Fig. ^ we show the ratio of the isovector ax- 
ial and vector charges as a function of the pion 
mass squared. The full-QCD results (diamonds) 
are consistent with the old quenched results (cir- 
cles) which are also shown: no mass dependence 
is seen. The weighted average of the three full- 
QCD points is g A /g v =1.21(8) and is consistent 
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Figure 4. The u (square) and d quark (cross) con- 
tributions to the d\ matrix element in full QCD. 

with the quenched results. 

In Fig. we present the quark momentum 
fraction in quenched (circle) and full QCD (di- 
amond). Though not yet renormalized, in the 
quenched approximation it does not show any 
curvature toward the chiral limit. The first mo- 
ment of the transversity is shown in Fig. [3J In 
both the quenched and full QCD we obtain the 
chiral limit of the bare matrix element by linear 
extrapolation. We use the RI-MOM scheme with 
one loop running to non-perturbatively renormal- 
ize the matrix element. The quenched result is 
1.194(25) (MS, 2 GeV) and the full QCD result 
will be reported soon. 

We have also calculated the c?i matrix element 
for which lack of chiral symmetry results in a 
mixing with a lower dimensional operator. We 
reported last year that this mixing is suppressed 
with domain wall fermions in the quenched ap- 
proximation [5]. Our full-QCD d\ results, shown 
in Fig. 0] are again in sharp contrast with the un- 
subtracted Wilson fermion results indicating that 
DWF suppress the power divergent mixing. 

To summarize, the quenched calculations with 
DBW2 and DWF works well for nucleon matrix 
elements. The axial to vector charge ratio, g A / g v , 
is well reproduced to within 5 % of the experi- 
mental value even on a modest volume and rel- 
atively heavy quark mass. The dependence on 
quark mass is very mild. Likewise the moments 
of structure functions (x) q , d\, and (l)sq all show 
mild linear dependence on quark mass, though 
some of them are yet to be renormalized. The 
full-QCD calculations of these quantities seem to 
work as well. Similar mild linear dependence on 
the quark mass was observed, extrapolating to es- 



timates not contradicting with the quenched ones 
(where renormalization is done), albeit with lim- 
ited statistics. With QCDOC computer JHj) we 
will be adding more observables such as flavor- 
singlet quantities and form factors and higher mo- 
ments of structure functions that require finite 
momentum, and moving to larger lattice volumes 
and smaller quark masses. 
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